Andreev bound states and tunneling characteristics of a 
non-centrosymmetric superconductor 
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The tunneling characteristics of planar junctions between a normal metal and a non-centro- 
symmetric superconductor like CePtsSi are examined. It is shown that the superconducting phase 
with mixed parity can give rise to characteristic zero-bias anomalies in certain junction directions. 
The origin of these zero-bias anomalies are Andreev bound states at the interface. The tunneling 
characteristics for different directions allow to test the structure of the parity-mixed pairing state. 
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Since the early sixties tunneling spectroscopy has 
played an important role in gathering information about 
the gap function of conventional superconductors [l[ . In 
the context of unconventional superconductivity tunnel- 
ing appeared as a tool to probe the internal phase struc- 
ture of the Cooper pair wave functions. Surface states 
with sub-gap energy, known as Andreev bound states, 
provide channels for resonant tunneling leading to so- 
called zero-bias anomalies. Zero-bias anomalies observed 
in high-temperature superconductors showed the pres- 
ence of zero-energy bound states at the surface, giving 
strong evidence for (i-wave pairing @, H, 0, @, @| . Simi- 
larly the tunneling spectrum observed in Sr 2 Ru04 is con- 
sistent with the existence of chiral surface states as ex- 
pected for a chiral p- wave superconductor 0,H,@, I3> 111- 
Thus quasiparticle tunneling has emerged as an impor- 
tant phase sensitive probe for unconventional supercon- 
ductors. 

The recent discovery of superconductivity in the heavy 
Fermion compound CePtsSi has motivated intense exper- 
imental and theoretical studies, since its crvstal symme- 
try is lacking a center of inversion 
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was shown that the presence of antisymmetric spin-orbit 
coupling in such a material could be responsible for a 
number of intriguing phenomena, such as the rather high 
upper critical field exceeding the standard paramagnetic 
limit 
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13|| . The absence of inversion symmetry yields 
a classification scheme of the pairing symmetry, different 
from the standard distinction between even- and odd- 
parity states. Actually, the pairing states in these non- 
centrosymmetric superconductors can be viewed as states 
of mixed parity which do not have a definite spin-singlet 
or spin-triplet configuration, either. For CePtsSi various 
experiments can be interpreted in a way that the pairing 
state has the highest possible symmetry, but would de- 
velop line nodes in the gap due to a mixing of the s-wave 
and a p-wave component. This state belongs to the A\- 
representation of the generating tetragonal point group 
C± v for CePt 3 Si. While thermodynamic quantities give 
good evidence for this behavior, it would be desirable 



to have additional proof for this kind of state through a 
phase sensitive test. The characteristics of quasiparticle 
tunneling may provide such a means. Despite the fact 
that the parity mixing Ai-phase has the full symmetry 
of the system, it gives rise to characteristic non-trivial 
features in the tunneling spectra (cf. also [17]). 

In the present study, we address the problem of 
possible Andreev bound states at the boundaries of a 
non-centrosymmetric superconductor like CePtaSi and 
demonstrate that quasiparticle tunneling could give im- 
portant information on the gap structure of such a ma- 
terial. The theoretical framework for our calculations is 
given by the quasiclassical Eilenberger theory [l8 , H 2(| . 
Our starting point is the quasiclassical bulk propagator 
gB in a non-centrosymmetric superconductor with an- 
tisymmetric Rashba-type spin-orbit coupling, which has 
already been derived in Ref. [lq : 



g B = -«7T 
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where the upper two components are given as 
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and similar relations hold for the lower components / 
and g. On the right hand side of Eqs. ([2]) and |[3]), u> n = 
(2n+l)7rfcsT denotes the Matsubara frequency, and both 
the gap amplitudes A± and the spin matrices a± are In- 
dependent. For a spherical Fermi surface parametrized 
in the standard way by k = (cos <f> sin 9, sin (f> sin 6, cos 9) , 
the gap amplitudes are 



A± = ip ± A sinff, 
while the spin matrices may be written as 



1 



a± = 2 V ±ie« 
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The variables if; an( A A are taken to be real and denote 
the gap amplitudes of the s- and p-wave components of 
the parity-mixed superconducting state. Altogether, this 
result is interpreted as a superconducting phase on two 
split spherical Fermi surfaces, each of them having a non- 
trivial gap function with amplitudes given by Eq. Q 
and complex spin structures according to Eq. ([5])- For 
technical simplicity we assume here that the splitting is 
small, such that both Fermi surfaces have identical size 
and a description by a one-band formalism is possible as 
described by Eqs. {U-©. A generalization to a really 
split Fermi surface is straightforward, but requires a more 
extended notation without giving additional insights. 

In order to finally include the boundary effects, it is 
convenient to use the Riccati-parametrization of the qua- 
siclassical theory [21, 22 1. Within this technique, the full 
quasiclassical propagator g in combined particle-hole and 
spin space may be written in terms of two coherence func- 
tions 7 and 7 in the following way [22l |: 
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The coherence functions 7,7 themselves are 2 x 2-spin 
matrices, which contain local information about the 
particle-hole coherence for a given k-direction. Partic- 
ularly, 7s and jb in the bulk of a non-centrosymmetric 
superconductor can be obtained just by comparing the 
known result of Eqs. fl])-© with the general form in Eq. 
©, leading to: 



7S = lil + g)- 1 ! 
1b = ij-\X-g), 



which can be written as 



IB = -(7+ct + +7_<3-_)<7 2/ 
IB = y {l+a+ +7_<7_), 



where we define 
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The advantage of this Riccati formalism is, that spatial 
inhomogeneities such as the boundaries can be taken into 
account in a straightforward way 2^, 23|, 24 1 . 

In the vicinity of the boundary the pair potential can 
develop a spatial dependence and deviate from the bulk 
form. Since we focus mainly on the qualitative aspects, 
however, we ignore this effect in the following and assume 
the gap is constant throughout the superconductor. The 
usual procedure within the Riccati formalism, namely 
solving differential equations of the Riccati type numer- 
ically for given k-direction, is unnecessary in this case. 
The presence of the surface enters through rather sim- 
ple boundary conditions for the coherence functions [22l | . 



For an impenetrable interface, the coherence functions 
7s j 7s at the surface of the superconductor are given an- 
alytically in terms of the bulk coherence functions: 

7s(k ou t) = 7s(k in ) 7s(k ou t) = 7s(k ou t)- (10) 

Here ki n and k out denote the incoming and outgoing di- 
rections of the wave vector, which are connected through 
the condition of specular scattering. Now the full quasi- 
classical propagator gs at the boundary is easily gen- 
erated from the coherence functions according to Eq. 
(|10j) . From gs,ii we obtain, for example, the angular 
resolved local density of states at the boundary of the 
non-centrosymmetric superconductor: 



N(E,k) 
N ^Re 



(11) 
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where E denotes the quasiparticle energy and Nq is the 
density of states of the normal state. With (...) denoting 
an average over the Fermi surface, the local density of 
states is just given as N(E) — (N(E,k)). Note, that 
at the boundary it is sufficient to average over half of 
the Fermi sphere only, i.e. we can restrict ourselves to 
quasiparticles traveling in k out directions. 

Eventually, we also obtain the normalized differen- 
tial tunneling conductance for a normal metal/non- 
centrosymmetric superconductor junction at low temper- 
atures T <§; X" c . To first order of the interface transmis- 
sion probability, it can be generated from the angular 
resolved local density of states of Eq. (fTTj) in the follow- 
ing way [iH| : 

Gr(eV) = (cosa D(a)N(eV,k))(cosa D(a)}- 1 . (12) 

The function D(a) denotes the transmission coefficient 
of the interface depending on the angle of incidence a, 
which is related to the normal vector n of the boundary 
by cos a = n • k out . In accordance with the result for 
an interface of the 5-function type 0, the transmission 
coefficient is taken to be 



D(a) 



D cos 2 a 
1 — Dq sin a 



(13) 



where the parameter Dq specifies the maximum transmis- 
sion coefficient for the angle of normal incidence a — 0. 

For the surface with normal vector n || z the rela- 
tion between k; n and k out needed to evaluate 75 (km) 
in Eq. (JTOj) is given by in = out and 6 m = it - 6> out , 
and the angle of incidence is a — out - In this situa- 
tion, we find immediately that A±(kj n ) = A±(k out ) and 
cr±(k in ) = <7±(k out ), which directly yields 7s(k ou t) = 
7s (kout)- Hence, both 7s and 7s are finally given by 
the bulk coherence functions for the same k-direction, 
and therefore the resulting angular resolved density of 
states N(E, k) is identical to the one of the bulk. No 
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FIG. 1: Normalized low temperature conductance G(eV) for a 
normal metal/non-centrosymmetric superconductor junction 
with boundary orientation n || z. The ratio q of the gap 
amplitudes is varied from to 1.5 in steps of 0.1. Between 
two subsequent curves, a vertical offset of 0.2 has been used 
for clarity of the plot. The transparency of the interface is set 
to D = 0.1. 



FIG. 2: (color online) Typical angular resolved density of 
states N(E, k) at the boundary of a non-centrosymmetric su- 
perconductor with normal vector n || y. For this example, we 
set q — 0.5 and 9 — it/2, i.e. we deal with k-vectors of the 
a;,j/-plane only. Bright colors represent high spectral weight. 
The branches of low-energy Andreev bound states according 
to Eq. (|14|l are clearly visible. 



sub-gap Andreev bound states are generated by surface 
scattering. For small but finite energies the shape of the 
local density of states is, however, determined by the ac- 
tual nodal topology of the gap functions given in Eq. 
This behavior is then passed on to the resulting tunneling 
conductances, which are presented in Fig. [T]for different 
values of the ratio q — if)/ A. Starting with q — 0, we have 
a pure p-wave state with point nodes on the poles of the 
Fermi surface. This leads to a parabolic behavior for low 
bias voltage reflecting the independence of the density 
of states. When q is increased, the gap function A + of Eq. 
(|4]) becomes fully gapped (albeit anisotropically) , while 
A_ exhibits immediately line nodes developing around 
the top and bottom of the Fermi sphere. Accordingly, 
a linear dependence can be observed for low bias. For 
q = 1 the line node of A_ has moved to the equator of 
the Fermi sphere. This removes the linear low-bias be- 
havior in G, since the nodal quasiparticles traveling par- 
allel to the surface cannot contribute to the conductance 
Eq. any longer. As soon as q > 1 however, also the 
gap function A_ is nodeless, which results in the open- 
ing of a corresponding gap in the conductance situated 
between the normalized energies Ej A = ±|g — 1|. 

Next we turn to the boundary with normal vector n _L 
z where we can choose freely n || y. Here, the relation be- 
tween incoming and outgoing quasiparticle directions is 
determined by (p- m = — </w and 6- m — 6> out . Furthermore, 
the angle of incidence is a = arccos(sin out sin# out ). In 



contrast to n || z, low-energy Andreev bound states can 
exist for this orientation. By examination of the angular 
resolved local density of states N(E, k) given in Eq. (jTTJ) 
we find the condition that these bound states only occur 
for k-vectors with A_ < 0, i.e. the contributions come 
from those parts of the Fermi surface, where the p-wave is 
stronger than the s-wave (sin# > q). Actually, for every 
angle 9 with A_ < there is a whole branch of bound 
states living between the energies ±|A_| determined by 



E 2 - A + A_ - 



COS( 



= ±- 



A\ -E 2 



A 2 _ - E 2 



(14) 



Zero-energy Andreev bound states occur for quasiparticle 
directions with <j) = n/2 which corresponds to vanishing 
transverse momentum along the x-direction. An example 
of the low-energy Andreev bound states in the plane 9 = 
7r/2 for q — 0.5 is provided by Fig. [2] 

The normalized low-temperature conductance for this 
situation is presented in Fig. [3l We find an extended 
zero-bias conductance peak for q = 0, analogous to the 
one for a pure p-wave superconductor with a cylindrical 
Fermi surface @, [|| . This peaked structure at low bias 
voltages persists as long as q < 1 but becomes progres- 
sively narrower for q — > 1 . These structures are easily un- 
derstood in the context of the Andreev bound states. For 
a given q < 1 there is always a region around the equator 
of the Fermi sphere with sinf? > q, where A_ < 0. Quasi- 
particles from those parts of the Fermi surface support 
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FIG. 3: Normalized low temperature conductance G(eV) for a 
normal metal/non-centrosymmetric superconductor junction 
with boundary orientation n _L z. The ratio q of the gap 
amplitudes is varied from to 1.5 in steps of 0.1. Between 
two subsequent curves, a vertical offset of 0.2 has been used 
for clarity of the plot. The transparency of the interface is set 
to D = 0.1. 

low-energy Andreev bound states up to normalized ener- 
gies of ±| sin — q\ according to Eq. (fl"4"|) . The maximum 
possible energy ±|1 — q\ originates from trajectories at 
the equator and corresponds to the width of the zero-bias 
anomalies shown in Fig. [3J Moreover, the total weight 
around zero bias naturally decreases when q is increased, 
because the region on the Fermi sphere with A_ < 
is reduced and accordingly the contributions to Andreev 
bound states. If q > 1, however, A_ > everywhere on 
the Fermi surface and Andreev bound states do not exist 
anymore. Therefore, the zero-bias anomaly disappears in 
the tunneling conductance. 

The results presented above have been achieved for a 
typical tunnel junction with small interface transparency 
Do = 0.1. Tunnel junctions with other interface trans- 
parencies chosen not too large to be consistent with the 
validity of Eq. (fl~2)) show the same characteristic behav- 
ior: For boundary orientations n _L z and q < 1 the tun- 
neling conductances exhibit a strong zero-bias anomaly. 

Our discussion shows that the tunneling characteristics 
could give important information on the gap structure 
of non-centrosymmetric superconductors, in particular in 
view of the parity-mixing. Important is also the direc- 
tion dependence which would provide an important test 
for the theoretical picture developed for CePtaSi sofar. It 



has also been anticipated that the ratio q could be tem- 
perature dependent. Thus it is possible that the tunnel- 
ing spectrum displays qualitative changes as a function 
of temperature. 
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